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The problem of calculating the four-nucleon bound state properties for the case of realistic two- 
and three-body nuclear potentials is studied using the hyperspherical harmonic (HH) approach. 
A careful analysis of the convergence of different classes of IfH functions has been performed. A 
restricted basis is chosen to allow for accurate estimates of the binding energy and other properties 
of the *He ground state. Results for various modern two-nucleon and two- plus three-nucleon 
interactions are presented. The *He asymptotic normalization constants for separation in 2+2 and 
1-1-3 clusters are also computed. 

I. INTRODUCTION 

Rapid progress has been made during the last few years in the quantitative study of the A — A nuclear systems. Ever 
increasing computer power, development of novel numerical methods, and significant refinements of well-established 
techniques have allowed the solution of the four-nucleon bound state problem with a control of the numerical error at 
the level of 10-20 keV (the experimental a-particle binding energy (BE) being 28.30 MeV), at least for Hamiltonians 
including only nucleon-nucleon (NN) interaction models [y. In the latter work, the BE and other properties of the 
a-particle were studied with the AV8' NN interaction, and the different techniques produced results in very close 
agreement with each other (at the level of less than 1%). 

In Refs. 1^1^, realistic potential models have been used to describe the a-particle bound state. Those potential 
models consist of the sum of a modern NN interaction plus a three-nucleon (3N) interaction. A modern NN interaction 
has the property of describing the NN database with a per datum close to one. Examples are the Argonne VIS 
(AV18) potential 0, the Nijmegen potentials ,6], the CD-Bonn potential 0,13 and the recently proposed potential, 
non-local in r-space, developed by Doleschall and coll. 0, |0 . This last potential has the remarkable property of 
reproducing simultaneously the NN bound and scattering data and the 3N binding energies. As is well known, the 
other models (AV18, Nijmegen and CD-Bonn) underbind the 3N system. Usually a 3N interaction is included in the 
Hamiltonian when these potentials are considered. The strength of the 3N interaction is properly tuned to reproduce 
the binding energy and this strength depends on the chosen NN potential. Examples of 3N interactions are the 
Urbana IX (UIX) 1\, the Tucson-Melbourne (TM) [H! and Brazil I?] potentials. From Ref. 3] we observe that all the 
NN-I-3N potential models which reproduce the deuteron and the 3N binding energies slightly overbind the a-particle. 
We further observe that the results obtained using different techniques jj, 0| for the AV18-f UIX potential model 
though close to each other are not in complete agreement. Clearly a clarification of these points would be welcome. 
Moreover, in recent years there has been a rapid pro gres s in developing new models of the nuclear interaction based 
on the application of the chiral perturbation theory [l^, [ij, [13 • From these studies one can hope to have a better 
understanding of the form of the NN and 3N interactions (the four-nucleon force is expected to be very small). All 
these potential models have to be studied in detail in the A — i and A = 4 systems. It is therefore very important to 
have powerful techniques for solving four-nucleon problems. 

The methods devised to tackle the problem of the solution of the non-relativistic Schrodinger equation 

i7*(l,2,3,4)-i?*(l,2,3,4) , (1.1) 

where H is the four-body nuclear Hamiltonian, are very different. In the Faddeev-Yakubovsky (FY) approach 
111 111 El. Eq. is transformed to a set of coupled equations for the FY amplitudes, which are then solved 

directly (in momentum or coordinate space) after a partial wave expansion. In the Green Function Monte Carlo 
(GFMC) method '?,'T^ one computes exp(-Ti7)$(l, 2, 3, 4), where $(1,2,3,4) is a trial wave function (WF), using 
a stochastic procedure to obtain, in the limit of large r, the exact ground state WF ^E*. These two techniques have 
also been applied to the case where the nuclear Hamiltonians includes a 3N interaction. The Stochastic Variational 
Method (SVM) jMIHI and the Coupled Rearrangement Channel Gaussian-Basis method (CRCG) jUjlll provide a 
variational solution of Eq. (|l.l|l by expanding the (radial part of the) WF in gaussians. The two techniques differ in 
the way they determine the non-linear coefficients of the expansion: in the SVM random choices are used to select the 
optimum set, whereas in the CRCG technique the non linear coefficients are chosen in geometrical progression in such 
a way that only a few of them have to be varied. Very recently two other new techniques have been proposed. In the 
no-core shell model (NCSM) method 0jll3 the calculations are performed using a (translational-invariant) harmonic- 
oscillator (HO) finite basis P and introducing an effective P-dependent Hamiltonian Hp to replace H in Eq. (|l.l|l . 
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The operator Hp is constructed so that the solution of the equation Hp^{P) = Ep'^{P) provides eigenvalues which 
quickly converge to the exact ones as P is enlarged. The effective interaction hyperspherical harmonic (EIHH) 
method jl^ is based on a similar idea, but the finite basis P is constructed in terms of the hyperspherical harmonic 
(HH) functions. 

In the present work we would like to address the problem of calculating the a-particle properties, using a nuclear 
Hamiltonian containing modern two- and three-nucleon interactions, by expanding the WF in terms of the HH 
functions. Our intention is to obtain converged binding energies at the level of 20-30 keV. The motivation is twofold. 
First we would like to reduce the theoretical error in the determination of the a-particle bound state properties. 
Second, the HH techniques can also be extended to treating four nucleons scattering states, as has been possible for 
the j4 = 3 system ^| using a similar technique. This program is currently under way and a preliminary report has 
been already published 2^ . The richness of phenomena in the four nucleon scattering and reactions will be an ideal 
laboratory for studyingand testing newer models of the nuclear interaction. 

In an earlier work j29| , the authors determined the solution of Eq. variationally by expanding the WF in a 

basis of correlated Hyperspherical Harmonic (CHH) WF's. The space part of such a basis consisted of products of 
correlation factors F and HH functions. The correlation factors F were chosen so as to take into account the strong 
correlations induced by the NN potential, especially at short inter-particle distances. The introduction of such factors 
substantially improved the convergence of the expansion. This made it possible to obtain reasonable estimates for 
the ground state energy of the a particle and some selected observables in n — ^H and p— '^He elastic scattering using 
a rather limited basis set |29, 30, 31] . However, due to the complexity of F, the spatial integrations were performed 
by using quasi-random number techniques. The precision of the required matrix elements was therefore limited, and 
the inclusion of a greater number of states was problematic. 

When the four-nucleon WF is expanded in terms of the uncorrelated HH basis (i.e. setting F = 1) most of the 
integrations can be performed analytically, and the remaining low-dimensional integrals can be evaluated by means of 
efficient quadrature methods. However, due to the particular structure of the NN potential, which is state dependent 
and strongly repulsive at short distances, a very large number of basis elements are required. For that reason, the 
application of the HH technique to studying the A — A nuclear system has encountered serious convergence problems. 
Few four-body HH calculations have been attempted so far for realistic interactions J32^ . Even for central 

or super-soft-core potentials the problem of the slow convergence of the HH expansion has not been completely 
overcome [s^, 13 US ■ The reason for these difhculties is related to the slow convergence of the basis with respect to 
the grand angular quantum number K and to the large number of HH states with a given K . For example, for an 
accurate description of the a-particle ground state, antisymmetric spin-isospin-HH states vlp to K = 60 have to be 
included. However, the number of such states already for K = 2Q is greater than 1,000 and it increases very rapidly 
with K. It is therefore clear that a hrute force application of the method is not possible even with sophisticated 
computational facilities. 

The approach analyzed in the past was to select a suitable subset of states [s^ 113, 113 . In those papers it resulted 
quite clearly that the quantum number K is not the unique parameter important for studying the convergence of the 
basis. Let us recall that a 4-body HH function is specified by three orbital angular momentum quantum numbers £i, 
£2, ^3 and two additional quantum numbers 712, (which are non-negative integers) related to the radial excitation of 
the system. The grand angular quantum number is defined tohe K = £1 -1-^2 +^3 +2(^2+^3). Note that C = (1+^2+^3 
and K are even (odd) numbers for positive (negative) parity states. In Ref. i36i | , the basis was restricted to including 
HH states with a few choices of ii, £2, is values, and large values of n2 and n^. The calculations performed 32„ .33| 



were however limited by the computer power available at that time. In this paper, it is shown that HH states having 
£ = £1 + £2 + < Q are sufficient in order to obtain a four digit convergence. However, the number of HH states with 
C< 6 is still huge and additional criteria for selecting a reduced basis have to be specified. 

It is possible to organize the HH states in terms of the number of particles correlated. For example, there is a class 
of elements which depends only on the coordinate of two particles, the so called Potential Basis (PB) [s^- Such a 
basis therefore takes the two-body correlations into account. However, even in the case of simple model interactions, 
the BE's B obtained by restricting the expansion basis to the PB were found to be rather far from the exact values. 
For example, for the Malfliet-Tjon V (MT-V) central potential 48j, B calculated with the PB is approximately 1 MeV 
smaller than the exact value. For a realistic potential the situation is noticeably worse. However, it is clear that the 
procedure of classifying the HH states in terms of the number of correlated particles can be useful for distinguishing 
the importance of the various expansion terms. 

In the present paper the application of the HH expansion basis is developed by taking advantage of both strategies 
discussed above. Namely, HH states of low values of £1, £2, £3 are included first. Among them, those correlating 
only a particle pair are included first, then those correlating three particles are added and so on. In practice, the 
HH states are first divided into classes depending on the value of £ and n2, na. Let us denote with Mi{Ki) the 
number of states belonging to a class i with K < Ki. Such a number of states rapidly increases with Ki, in general 
A4i{Ki) « Mi{KiY for large Ki, where the numbers Mi are a set of constants and r = 1 or 2. The first and most 
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important classes should contain only a small subset of HH states ( "small" classes) , namely their Mi should be small, 
let us say Mi \. It is then relatively easy to include HH states of large K belonging to these classes. The classes 
containing successively larger numbers of states {Mi 3> 1) should be chosen possibly so as to contribute lesser and 
lesser to the expansion in such a way that the corresponding expansion can be truncated to small values of K. A 
careful analysis of the convergence properties of the various HH components has allowed for an optimal choice of the 
classes, so that accurate calculations of the a-particle properties could be achieved. 

An important aspect of a successful application of the HH method is related to the computation of the coefficients 
for the transformation of a HH function corresponding to a generic permutation of the 4 particles in terms of those 
constructed for a given permutation. Various approaches have been devised to deal with this problem 0, |^ [i^, 
0,0,^3- The usefulness of these coefficients is twofold. First, it is easy to identify the linearly dependent states 
and to avoid their inclusion in the expansion basis. The removal of these "spurious" states, which disappear after a 
proper antisymmetrization of the basis, is very useful as the number of linear independent states is noticeably smaller 
than the full degeneracy of the basis. Secondly, the matrix elements of a local two-body (three-body) potential energy 
operator are easily reduced to one-dimensional (tri-dimensional) integrations, which can also be performed beforehand 
and stored on computer disks. The matrix elements of non-local operators can also be reduced to low-dimensional 
integrals. The kinetic energy operator is easily obtained analytically. 

This study is the continuation of the application of the HH expansion to the three nucleon system performed in 
Ref. . Another possible extension of these studies are related to the application of the HH technique to heavier 
systems. In particular, we can point out that also for A> A the calculation of the multidimensional integrals related 
to the matrix elements of a local NN (3N) interaction reduces to a one-dimensional (three-dimensional) integration. 
The only difficulty in extending the method to heavier systems is the choice of a suitable and optimized subset of 
HH functions. We hope that the criteria used here to select an optimal subset of the basis could also be applied for 
systems with A> A. 

This paper is organized as follows. In the next section, a brief description of the properties of the HH functions is 
reported. In section 3, the choice of the basis is presented. The results obtained for the BE and other properties of 
the a-particle are presented in Sect. 4. Finally, the last section is devoted to the conclusions and the perspectives of 
the present approach. 



II. THE HH EXPANSION 

For four equal mass particles, a suitable choice of the Jacobi vectors is 

^ n + "^3 + rk 



ip ~ \l 2 ^'"'"^ 3 ^ ' 

I., = i'.irk-^^), (2.1) 



where p specifies a given permutation corresponding to the order i, j, k and m of the particles. By definition, the 
permutation p = 1 is chosen to correspond to the order 1, 2, 3 and 4. 

For a given choice of the Jacobi vectors, the hypcrspherical coordinates are given by the hyperradius p, defined by 



Cip + ?2p + Cip ' (independent on p) , (2.2) 

and by a set of variables which in the Zernike and Brinkman |38ll47| representation are the polar angles = {9ip, (pip) 
of each Jacobi vector, and the two additional "hypcrspherical" angles ip2p and ip^p defined as 

cos ip2p ^ ^ ^"^P cos yj3p = = ^ — , (2.3) 

Cip + i'ip y Cip + ^2p + Cip ^ 

where £,jp is the magnitude of the Jacobi vector ^^p. The set of the variables ^ip,£,2p,£,3p, 'P2p, 'Psp is denoted hereafter 
as J7p. To simplify the notation for p = 1, the subscript "1" will be sometime omitted. The expression of a generic 
HH function is 



yei',e2,e3,L2,n2,n3i^p) 



(y^Alp)yi2{^2p)) ^ Y^.ihp) 



V'„\f^fHv2p,^3p) . (2.4) 



"2,13 

LM 
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where 



"2, "3 



sin^^ ip2p cos^^ ip2p sin^^+^^+^"^ ip^p cos^ 



f3p 



^ (cos 2.^2p)i^3 (C0s2(^3p) , 



and P"'** are Jacobi polynomials. The coefficients M^^'^^f^ are normalization factors, given explicitly by 

3 r 



■'^n2,ii3 



n 

J=2 



r(^.-«.-^.-V2)r(n, + ^, + 3/2) 



where Vj = iiTj + (3j — 5)/2 with if^ defined to be 



if2 = ^1 + ^2 + 2n2 , ^3 = + ^3 + 2n3 = X 



(2.5) 



(2.6) 



(2.7) 



and K is the grand angular quantum number. 

The HH functions are eigenfunctions of the hyperangular part of the kinetic energy operator A^. In fact, for j4 = 4 
the latter operator can be written using the variables {p, fip} as follows 



i=i,3 



8d_ 
pdp 



and 



A2(f7p) + K{K + 7))3^,^.:tl,L„„„„3(f^p) = 



(2.8) 



(2.9) 



Another important property of the HH functions is that /o^iV^'^'^ „^ n^^^p) homogeneous polynomials of the 
particle coordinates of degree K. 

The WF of a state with total angular momentum J, parity tt and total isospin T can be expanded over the following 
complete basis of antisymmetrical hyperangular-spin-isospin states, defined as 



KLSTJtt 



12 



KLSTJtt 



{i,j;k;Tn) 



where the sum is over the 12 even permutations p, and 



KLSTJtt 



jKLM 



{ij]k;m) 



yK.LM :^ 



SJ JJ^ 



(2.10) 



(2.11) 



TT^ 



Here, 3^/^ £^'£3 L2 "2 n3(^p) ^^e HH state defined in Eq. (|2.4|l . and (t.^) denotes the spin (isospin) function of 
particle i. The total orbital angular momentum L of the HH function is coupled to the total spin S to give a total 
angular momentum J, J^. The quantum number T specifies the total isospin, while tt = (— l)^i+^2+<'3 jg ^^j^g parity of 
the state. The integer index labels the possible choices of hyperangular, spin and isospin quantum numbers, namely 



A* = {h,<^2,<^3,L2,n2,n3, Sa, Sb,Ta,Tb} , 



(2.12) 



compatible with the given values of K, L, S, T, J and tt. Another important classification of the states is to group 
them into "channels": states belonging to the same channel have the same values of angular £1,12,^3, L2, L, spin 
Sa, Sb, S and isospin Ta,Tb,T quantum numbers but different values of rt2, n.3. 

Each state x^/KLSTJtt g^tering the expansion of the four-nucleon WF must to be antisymmetric under the exchange 
of any pair of particles. Consequently, it is necessary to consider states such that 



KLSTJtt 



{i,j;k;m) 



KLSTJtt 



{j,i;k;Tn) 



(2.13) 



Under the exchange i ^ j, the Jacobi vector changes its sign, whereas and remain unchanged, and, 
therefore, the HH function y^^^f^\^ 713 (^p) transforms into itself times a factor (—1)^^ (see Eqs. (|2.1f) and (|2.4|) '). 
On the other hand, the spin-isospin part transforms into itself times a factor (— 1)-S£i+Ta j-qj. ^j^g i j exchange. Thus, 
the condition H2.13|l is fulfilled when 

^3 + 5„ + r„ =odd . (2.14) 
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The number Mklsjt-k of the antisymmetrical functions vJ/^-^^S-^Tt having given L, S, T, J and tt values but 
different combination of the quantum numbers fi is in general very large. In addition to the degeneracy N^l-k of the 
HH basis, the four spins (isospins) can be coupled in different ways to S (T). However, many of the states \ij^^SJTTr 
are linearly dependent amongst themselves. In the expansion of a four-nucleon WF it is necessary to include the subset 
of linearly independent states only. To search for the independent states, the essential ingredient is the knowledge of 

^^mI^'^"" = (^^<^STJ.^^KLSTJ.^^ ^ (2.15) 

where ()n denotes the evaluation of the spin-isospin traces and the integration over the hyperspherical variables. 

The calculation of the matrix elements of the Hamiltonian is considerably simplified by using the following trans- 
formation 

$^^^"^(*,j;A;m) = 5:a™-(p)$^^^"-(l,2;3;4) . (2.16) 

The coefhcients o.^^^'^''^ {p) have been obtained using the techniques described in Ref. The states \^^^STJit 
be written as 

^KLSTJ. ^ ^KLSTJ.^KLSTJ.f^^^ 2; 3; 4) , (2.17) 

where 

^^^""^=E<^'"'^(^')- (2-18) 



The matrix elements of the norm can be easily obtained using the orthonormalization of the HH basis with the result 
that: 

^^.1.' - T.(<'J^'^yKu'^'^ ■ (2-19) 

fj." 

Clearly, 

^i^KLSTj^^^K'L's'T'j'^'^^^ = , if {KLSTJt:} ^ {K'L'S'T'J'tt'} . (2.20) 

Once the quantities N^fjf^'^'^^ are calculated, the Gram-Schmidt procedure can be used, for example, to eliminate 
the linear dependent states between the various •^^'^stj-k functions. 

We have found that the number of independent states M'^^grpj^ for given if, L, S, T, J and tt is noticeably smaller 
than the corresponding value of Mklstj-k ■ To give an example, we have reported in Table^a few values of Mklstj-k 
and M'j^^grpj^ for the case J = 0, T = 0, 7r = + corresponding to the ground state of the a-particle. As can be seen 
from the table, the values of Mklloo+ are very large also for moderate values of K, but M^^^qq^ are usually much 
smaller. 

The total WF can finally be written as 

Vff ^ ^^^^I^^^M^KLSTJ. ^ (2.21) 

KLST fJ. P 

where the sum is restricted only to the linearly independent states. The expansion coefficients, which depend on the 
hyperradius, are determined by the Raylcigh-Ritz variational principle. By applying this principle, a set of second 
order differential equations for the functions u{p) a.re obtained. These equations and the procedure adopted to solve 
them has been outlined in the appendix of Ref. |46l |. In this way, a large number of equation can be solved. 

The main problem is the computation of the matrix elements of the Hamiltonian. The kinetic energy operator 
matrix elements are readily calculated analytically, whereas the matrix elements of a local NN potential can be 
obtained by one dimensional integrations. To this aim, it is convenient to write the basis in the jj coupling scheme 

^KLSTJ. ^ ^ j^KLSTJ. ^KTJ.^^^ ^. 3. 4) ^ (2.22) 
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where 



1,2; 3; 4) 



(i£3(6)(5lS2)s„)^.^(l£,(6)s3 



J2, 



J2 



(2.23) 



and B^^^^"^^ are related to the coefficients A^^^'^'^'^ via Wigner 3j and 6j coefficients. Now, the integer index 
labels all possible choices of 



= {n3j3,Sa,j3,n2,h,j2,J2,h,jl,Ta,Tb} 



(2.24) 



compatible with the given values of T, J and tt. 

In terms of the states S^^'^'^(l, 2; 3; 4), it is easy to compute the matrix elements of an NN potential. For example, 
the matrix element of the isospin-conserving part Vjci^, 2) of the NN potential 



(S™(1, 2; 3; 4)|V>c(l, 2)|S^;'^-^-(l, 2; 3; 4))o = 



K'T'J-nl 



unless 



If Eq. 



{ J3 , n2 , ^2 , J2 , ^2 , 4 , Ji , T,, Tb, T} = {j^ , n'^ , 4 , j'^ , , £'i , j J , , T{, , T'} 



(2.25) 
(2.26) 



is verified, then 

(S™(1, 2; 3; 4)|t/,c(l, 2)|S5'^'^-(1, 2; 3; 4))^ 



A/'»i;l^/^AC:lf / d^3(cos^3)^+^-^+^Ksin^3) 

73 / \ nfl+^2+2n2+2/3 + ^/ o \ 

^ S £' S' (PCOS^Ja) P„3 - (C0s2(y93) X 



N5+2£i+2£2+4n2 



XP^^ 



+£2+2n2+2,4 + ^ 



(cos 2(^3) , 



(2.27) 



where g ^, g,{r) is the isospin-conserving part of the NN potential acting between two-body states and 

2S The one-dimensional integral given in Eq. H2.26|l can be computed numerically with high accuracy. The 

case of the isospin-breaking part of the NN interaction is a generalization of the previous case: now we can have 
{r„rb,T}^^K,r,',T'} aswell. 

The 3N interaction matrix elements are more difficult to compute and the adopted procedure is detailed in the 
Appendix. 



III. CHOICE OF THE BASIS 



The main difficulty of applying the HH technique is the selection of a restricted and effective subset of basis states. 
In fact, although the number of independent states proves to be much smaller than the degeneracy Mklstj-k of the 
basis, the brute force application of the method, i.e., the inclusion of all HH states having K < Km in the expansion 
and then increasing Km until convergence, would be destined to fail. In fact, due to the strong correlations induced 
by the NN potential. Km ~ 60 are necessary in order to obtain a good convergence. However, already for values of 
K > 20 it is very difficult to find the linearly independent states via the Gram-Schmidt procedure due to the loss of 
precision in the orthogonalization procedure. 

However, it is possible to separate the HH functions into classes having particular properties and advantageously 
take into account the fact that the convergence rates of the various classes are rather different. As discussed in the 
Introduction, we expect that the contribution of the HH functions describing the two-body correlations to be very 
important ^8 ]. Another criterion adopted is first to consider the HH functions with low values of ii. 

An important quantity in the choice of the classes is Mi{Ki), namely the number of linearly independent antisym- 
metrical spin-isospin-HH states x^i^^STJtt belonging to a class i and having K < Ki. Only even parity states have 
been included in the construction of the a-particle WF, and thus the discussion hereafter will be limited to consider 
only even values for Ki. In general, for a class i, the value A4i{Ki) is zero for Ki < Kf , due to the fact that the 
linearly-dependent states have been removed from the expansion. For example, as should be clear by inspection of 
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Table HJ there is only one linearly- independent state \^klsoo+ ^h]^ K = 0. If this state is included in the first class, 
the other classes must have at least K"' = 2, etc. For Ki ^ K^, Aii{Ki) reaches a sort of "asymptotic" value, given 
by A4i{Ki) « Mi{KiY\ The choice of the classes has clearly to be optimized so that the convergence for the classes 
with large values of Mi and could be reached for relatively low values of K. The specific values of Mi and are 
discussed below. 

To study the a-particle ground state we have found it very convenient to choose as follows (T is the total isospin): 

1. Class CI. In this class the T = HH states belonging to the PB are included . For A — 4, the PB includes 
states of the first three channels reported in Table lU (the only channels with £i ~ £2 ~ 0) with 712 = 0. As can 
be seen from Eq. H2.4|l . the corresponding states depends only on ^3^, and cosip^p — S.3p/p = rij/p, and therefore 
contain only two-body correlations. For this class, Kf = 0. For Ki > 4, Aii{Ki) = (3/2)Ki. Then, this is 
a "small" class. As will be shown in the next Section, this is also the most slowly convergent class, but since 
Ml = (3/2) and ri = 1, it is not difficult to reach the desired degree of accuracy. 

2. Class C2. This class includes the T = states belonging to the same three channels as those of class CI, 
but with n2 > 0. These states therefore include also part of the three-body correlations. The first linearly- 
independent states of this class appear for K = 4, therefore K2 = 4. Moreover, A^2(^2) = (3/4)(fsr2)^ + C>{K2) 
for K2 ^ 1. This can be considered a "small" class, too, and states up to K2 = 40 have been included in the 
present calculation without difficulty. 

3. Class C3. This class includes the remaining T = states of the channels having £1+^2 + ^3 = 2. The 
corresponding 20 possible channels are reported in Table UTI in rows 4 — 23. In this case — 2 and /^^{Ks) — 
5(if3)^ + ©(/-Ca) for K3 » 1. This is a fairly "large" class, but with the necessary care states with K3 « 34 can 
be still included in the expansion. 

4. Class C4. This class includes T = states belonging to the channels with £1+^2+^3 = 4. There are 57 channels 
of this kind. In this case K2 = 8 and it follows that Mi{Ki) = {57/'i){K4)^ + 0{Ka) for > 1. This is a 
"large" class, but its contribution to the a-particlc BE is, though still sizable, far less important than the first 
three classes. States of up to k, 28 have been considered. 

5. Class C5. This class includes T = states belonging to the channels with £1 + i?2 + ^3 = 6. There are 109 
channels of this kind. In this case = 12 and for ifg > 1 we have M^{K^) = {im/A){K^)'^ + 0{Kz). This 
is a very "large" class, but it contributes very little to the a-particle BE, as we shall see. Therefore, we can 
truncate the expansion already at « 20. 

6. Class C6. This class includes the states having T > 0. We have included in the expansion all the channels 
of this kind with £1 -|- £2 + ^3 < 2 (45 channels). In this case K^^ — and for Kq ^ 1 we have Aif,{K(,) ~ 
(45/4) (i^e)^ + C>{Kq). Also the contribution to this class is very tiny. 

The states belonging to the classes C2 and C3 describe the most important three-body contributions to the WF. 
The classes C4 and C5 take into account the remaining three and four body correlations ordered with increasing 
values of £i + £2 + £3.- 

The convergence is studied as follows. First, only the states of class CI with K < Ki are included in the expansion 
and the convergence of the BE is studied as the value of Ki is increased. Once a satisfactory value of Ki = Kim is 
reached, the states of the second class with K < K2 are added to the expansion, keeping all the states of the class 
CI with K < KiM- Then K2 is increased up to K2M in order to reach the desired convergence for the BE. With 
some extra work, it is possible at this point to optimize the basis by removing some of the K < K2M states of class 
C2 which give very tiny contributions to the BE. The procedure outlined is then repeated for each new class. Our 
complete calculation includes about 8, 000 states. 

It should be noticed that in the present calculation only HH functions constructed in terms of the Jacobi vectors 
given in Eq. (|2.1|l . referred to as the set A, have been considered. As is well known, there is another possible choice, 
namely 

^2p \l^irk + - ri - rj) , (3.1) 

^3p = '"j ~ '''i ' 

hereafter referred as the set B of Jacobi vectors. Considering, for example, the a-particle ground state, the HH 
functions J^sot A of the set A are more appropriate for describing those contributions to the WF corresponding to 
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{3 + 1} clustering structures, namely ^He + n or +p. The HH functions J^set b constructed with the set B, should 
be more suitable for describing the {2 + 2} clustering structures, such as the d + d configurations. It is rather obvious 
that the inclusion of HH functions of both sets should speed up the convergence in constructing the full state of the 
system |2^. If the expansion of the WF is done over only a particular set, those configurations in which other 
clustering structures are important would be generally described with difficulty and a slow convergence would result. 

In the present calculation we have included HH states 3^sct A only, i.e. constructed with the Jacobi vectors of the 
set A, since this has been found to be sufficient to reach the desired degree of convergence. In fact, the full basis 
considered (classes C1-C6) is large enough to include all the possible independent states for K < 20. Additional 
linearly independent states constructed with the set B would appear only for K > 22. As will appear clear below, 
the contribution of states with > 22 not belonging to classes C1-C3 is rather small. Therefore, in the present 
calculation it is not necessary to introduce states of the set B. However in the present formalism there would be no 
particular difficulty in also including states constructed with the set B. 

IV. RESULTS FOR THE a PARTICLE GROUND STATE 

In this section, the results obtained for the ground state of the a particle are presented. The convergence of the 
HH expansion in terms of the various classes is examined in Sect. II V Al The results obtained for the BE and other 
properties for a number of different interaction models are reported in Sect. lIV'Bl The origin of the T > components 
in the a-particle ground state is discussed in Sect. II VO The effect of the truncation of the NN and 3N interactions 
is studied in Sect. HVDl Finally, the calculation of the various "^He asymptotic normalization constants is considered 
in Sect. lTVEl 

A. Convergence 

In order to study the convergence, we have considered three different interaction models frequently used in literature. 
The first calculation has been performed using the MT-V potential a central spin-independent interaction. The 
parameters defining this potential can be found in Table I of Ref. "SS] and we have used h? /m = 41.47 MeV fm^. This 
potential has been used for a number of benchmarks. It does not contain any non central components, but it retains a 
rather strong repulsion at short interparticle distance going like l/r. It is therefore rather challenging for a technique 
where the correlations are not built in. In the second example, we have considered the AVI 8 potential model '5^| which 
represent a NN interaction in its full richness, with short-range repulsion, tensor and other non-central components, 
charge symmetry breaking terms, and Coulomb and other electromagnetic interactions. In the third case, we have 
added to the AV18 potential the Urbana IX model of 3N interaction (AV18-fUIX model). For the latter two 
models we have used ft^/m = 41.47108 MeV fm^. 

We study the convergence as explained in the previous Section, and the results presented in table IIIII are arranged 
accordingly. For example, the BE B reported in a row with a given set of values of JsTi , . . . , Kq has been obtained by 
including in the expansion all the HH functions of class Cj with K < Ki, i — 1, . . . 

For the MT-V potential, we observe a slow convergence of the classes CI and C2 and fairly large values of K have 
to be used. On the other side, they give 96% of the total BE. The contributions of the other classes are extremely 
small. The class C3 increases the BE by an additional 0.08 MeV, and the class C4 by less than 0.01 MeV. Class C5 
gives a negligible contribution, and class C6 has not been included in the expansion since for this potential isospin is 
a good quantum number and there is not any mixing with T > components. The final value B ~ 31.347 MeV is in 
good agreement with the results found in literature, whose "average" has been reported in the last row of Table ITTll 
There is approximately 10 keV of missing energy due to the truncation of our expansion as will be discussed at the 
end of this Subsection. 

For the AV18 potential, the first two classes give important contributions but a large amount of BE is still missing. 
The inclusion of the third class increases the BE by more than 3 MeV but 0.8 MeV are still missing. Since the second 
and third classes take into account a large part of the contributions of the three body correlations, this means that 
also the four body correlation are important. These are related to the configurations where the clusterization 2 -|- 2 is 
important. In our calculation, such configurations are included when the classes C4 and C5 are taken into account. 
The number of the states of class C4 increases very rapidly with K4 but fortunately the convergence is reached around 
K = 24. The gained BE is almost 0.8 MeV. There are no linearly independent states of class C5 with AT < 14 and 
its contribution is rather small. The convergence is again obtained around K = 24, but the gain in energy is only 
about 0.02 MeV. Since the number of states of this class is very large, for example A4q{20) w 800 when confronted 
with a very tiny gain in BE, a selection of the states has to be performed to save computing time and to avoid loss of 
numerical precision. For example, all the channels of class C5 with a total orbital angular momentum L — have not 
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been included in the expansion since their contribution is absolutely negligible. With their inclusion the procedure of 
Ref. 46J for finding the eigenvalue would become numerically instable. 

From table IIIll one can try to estimate the contribution of the states with £ = + £2 + ^3 = 8. From the previous 
discussion we have already seen that the states having C — A (class C4) contribute by about 0.8 MeV, while the states 
with £ = 6 (class C5) contribute by less than 0.03 MeV. Therefore, the states with £ — 8 are expected to give a 
negligible contribution to the a-particle BE. Finally, the inclusion of states with T > (class C6) increases the BE 
by another 20 keV, approximately. 

The convergence rate when considering the UIX 3N interaction is similar to the AV18 case. The corresponding 
results are reported in the last column of table IIIII (they have been obtained in the approximation described in 
Subsec. flV mi . Since the models most frequently used for the 3N interactions are rather soft at short interparticle 
distances, the convergence rate of the CI and C2 classes does not change appreciably. However, the 3N potential 
has a very strong state dependence and the convergence of the C3-C5 classes are now slightly slower. For example, 
the gain in BE of the C4 class is about 0.8 MeV without any 3N interactions, and it becomes about 1 MeV when 
including the 3N interaction. Our final results for the AV18 and AV18+UIX models agree well with the FY results of 
Ref. "FTl reported in the last row of Table Hill The convergence properties for other NN and NN+3N potential models 
has been found rather similar to those showed in Table Hill 

Finally, let us comment about the convergence rate of the expansion as a function of the maximal grand ang ular 
quantum numbers Ki of the various classes of HH states included in our expansion. Previous studies |3^ l38L l39l 
have shown that the trend of convergence toward the exact BE depends primarily on the form of the potential. In 
particular, for potentials which are given as functions of rfj (as, for example, those given as a sum of gaussians) the 
increase of BE with Ki diminishes exponentially. On the other hand, for potentials given as a function of r^^ (as a 
sum of exponentials or Yukawians), the increase of BE decreases as where p is a positive integer number. 

The value of p is smaller for potentials of Yukawa type due to the 1 /r divergence at the origin, but may depends also 
on the class of the HH functions whose convergence is studied. It is important to determine the value of Km at which 
the convergence starts to behave as stated previously. The asymptotic behavior of the convergence should be reached 
for HH functions whose kinetic energy c>c {H^ /m)K{K + 7)/pn is much greater than the BE, where po is a value of 
the hyperradius p for which ^'(po) can be regarded as small |37l| . In our studies, we have found that the asymptotic 
falling begins for K„i « 30 40. 

In order to study the convergence behavior we have indicated with B{Ki, K2, K^, K4, K^,Kq) the BE obtained by 
including in the expansion all the HH states of the class CI with K < Ki, all the HH states of the class C2 having 
K < K2, etc. Let us compute 

Ai{K) = B{K, 0, 0, 0, 0, 0) - B{K - 2, 0, 0, 0, 0, 0) , (4.1) 
A2(i^) - B(ifiM,i^,0,0,0,0)-B(i^iM,if-2,0,0,0,0) , i^m = 72 , (4.2) 
A3(i^) = B(/fiM, K2M, K, 0, 0, 0) - B{KiM, K^m, K-2,0, 0, 0) , i^iM - 72 , K2M - 40 , (4.3) 

and so on. The values obtained for A^, j = 1, 3 are shown in Fig. ^ for the MT-V potential model, together with the 
curves {1/K)p for p = 5 (the curves have been constrained to fit the high K part of the Ai^3{K) values ). As can be 
seen in Fig.^ all the energy differences Ai, A2 and A3 decrease as 1/K^ for K > 20, approximately. However, for a 
given K, there is a clear hierarchy Ai{K) ^ A2{K) ^ A3{K). Note that there are slight fluctuations in the A{K) 
as K is increased (this is evident in particular for A3). 

The values obtained for A^, i = 1,4 for the AV18 potential are reported in Fig. 13 The decrease of Ai, A3 and 
A4 clearly follows a law {1/K)p with p = 7. The behavior of the energy difference A2 can be approximated either 
by a 1/K^ or a 1/K'^ law. The faster decrease of these results compared to the previous case is due to the fact that 
the AV18 potential does not diverge at the origin, while the MT-V has a 1/r divergence. The study of Ref. [S^l, in 
fact, predicts a difference of two units in the exponential coefficient for the two cases (Yukawian potentials vs. regular 
potentials). For fixed K, also for the AV18 we note a systematic hierarchy Ai{K) A2{K) ^ A3{K) ^ A4{K), 
although less pronounced than in the MT-V case. The same behaviour is observed when the UIX 3N potential is 
included. 

From the observed simple behavior, we can readily estimate the missing BE due to the truncation of the expansion 
to finite values oi K — K. Let us suppose that the states of class i up to K = K have been included and to have 
computed Ai{K). Then, the missing BE due to the states with K = K + 2, K + A, . . ., is given by 

{AB),^c{K,p)A,{K) , c{K,p)= f-j , (4.4) 

K='K+2,'K+4:,... 

where c{K,p) is a numerical coefficient. For example, let us consider the "missing" energy for the Class CI in the 
MT-V case. In this case 7? = 72 and p = 5 and c(72,5) = 8.51. Since Ai(K = 72) = 0.99 keV, we find that 
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ABi = 8.4 keV. Adding this value to B{72, 0, 0, 0, 0, 0) we can extrapolate the BE for the case of the inclusion of the 
whole class CI: B{72, 0, 0, 0, 0, 0) + {AB)i « 30.041 MeV. The BE obtained corresponding to this case (converged PB 
expansion) has been computed very precisely in Ref. [H3| using a "pair-correlated" potential basis (PPB). In such an 
expansion, each PB function is multiplied by a pair correlation factor and this allows for a very rapid convergence of 
the expansion. In that paper, we found i?(PPB) — 30.042 MeV which agrees very well with the above extrapolated 
value. To reach such a value, it would be necessary to use K w 150 for the class CI. 

For the AV18, we find c(72, 7) = 5.52, Ai(K = 72) = 0.24 keV and therefore (AB)i w 1.3 keV, a rather tiny 
quantity. Due to the faster convergence for this potential like l/K"^, it does not seem necessary to increase Ki any 
further in this case. 

The "missing" energy of the other classes can be estimated in the same way. However, to estimate the "missing" 
energy for the whole calculation due to the truncation of the expansion of the first class up to if < Ki , of the second 
class up to if < K2, etc., we cannot simply add the {AB)i, z = 1, . . . , 6 so obtained. The reason is that, for example, 
the inclusion of the HH states of classes C2, C3, . . ., also alters the convergence of class CI, etc. by a small amount. 
To study the "full" rate of convergence, let us consider 

Ai(X) = B{K,K2M,K3M:K4M,K5M,KeM)- 

B{K — 2, if2M, K^M, KiM, K^Mt Kqm) , 

A2{K) = B(ifiM,if,if3M,if4M,if5M,if6M) - 

B{KiM,K - 2,K'iM,K4M,K^M,KQM) , (4.5) 
AaiK) = B{KiM,K2M,K,KiM,K5M,KeM)- 

B{KiM, K2M, K ~2, K4M, K-^M, Kqm) , 

and so on. Clearly Aq{K) = Aq{K). The differences between Ai{K) and Ai{K) ioi- i — 2-^ 5 have been found to be 
negligible. Only the differences between Ai{K) and Ai{K) are sizable. In any case the behavior of Ai{K) for large 
K is the same as that discussed for A{K). Therefore, we propose to estimate the "total missing" BE by using the 
formula 

(AB)t = J2 <K,M,P)A,{K^M) , (4.6) 

1=1,6 

where p — 5 {7) for the MT-V potential (AVIS and AV18+UIX). To give an example, the values for Ai{KiM) and 
c{KiM,p) computed for the MT-V case are reported in Table Hvl from which it is possible to derive that {AB)t ~ 11 
keV. If this value is added to 5(72,40,34,28,0,0) = 31.347 MeV, we obtain 31.358 MeV, which is in very good 
agreement with the results obtained by other groups. For the AV18 potential, Eq. 1)4. 6|l gives {AB)t = 12 kcV 
and if this value is added to 5(72,40, 34,28,24,16) = 24.210 MeV, we obtain 24.222 in close agreement with the 
FY estimates_of 24.25 MeV of Ref. (51| and 24.223 of Ref. |53|. Note that for the class C2 we have computed the 
coefficient c{K,p) with p = 7. Using p = 6, we have c(40,6) = 3.52 and (Ai?)2 — 3.20 keV (instead of 2.60 keV), a 
very small change. For the AV18-I-UIX model, the same procedure allows for an extrapolated BE estimate of 28.474 
MeV, again in agreement with the FY value 28.50 MeV. Note that the FY BE results are quoted with an uncertainty 
of 50 keV due to the truncated model space in their calculations jsj . 

B. Results 

The values obtained for a number of different potential models after including the states of the 6 different classes up 
to the values Ki = 72, K2 = 40, K3 = 34, K4 = 28, = 24 and Kq — 16 (the last two values only for the realistic 
cases, for central potential we have taken — Kq = 0) are presented in Tables El and IVII Table presents the 
results for some central potential models, while Table IVII reports the results for various realistic potentials with and 
without including different models for the 3N forces, too. The BE's obtained by using the extrapolation technique 
described in the previous section are enclosed in parentheses. Results obtained by other techniques are also reported. 

Let us consider first the central potentials (for all of them we have taken f? jm = 41.47 MeV fm^). We have 
selected 5 different potential models, i.e. the Volkov I55j| . the Afnan-Tang S3 (ATS3) [S^, the Minnesota 0|, the 
MT-V and the Malfliet-Tjon version I/III (MT-I/III) [ig. These potentials have been used by several groups to 
produce benchmark calculations, but unfortunately for some of them different versions exist. The parameters of the 
first 4 potentials mentioned above can be found in Table I of Ref. while the version of the MT-I/III used has the 
same parameters as reported in Table I of Ref. The Volkov and MT-V are spin-independent, while the other 3 
potentials are spin dependent. Note that it is customary to include the point-Coulomb interaction (e^ = 1.44 MeV fm) 
with the Minnesota potential, while the MT-I/III version acts only on s- waves. Clearly for this group of potentials, 
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the total orbital angular momentum is a good quantum number and therefore we have included in the WF's only the 
channels with L — 0. 

The first example is the Volkov potential with Majorana parameter M = 0.6. As can be seen in Table El our result 
agrees very well with the estimates by other techniques, especially with the one using the SVM 20] . The Volkov 
potential, given as a sum of gaussians, has a very soft core and therefore the induced two-body correlations in the 
ground state WF are weaker than in the other cases. In fact, we have found that the convergence of the HH expansion 
is in this case much faster (it is reached for Ki « 30). Since inclusion of HH states with fairly low values of the 
grand angular quantum number are sufficient to obtain convergence, a successful HH calculation for this potential 
was already possible 20 years ago [ssf 

Others central potentials often used in the literature are the ATS3 and Minnesota potentials. Both are given as a 
sum of gaussians but have a rather strong repulsion at short interparticle distances. This induces important two-body 
correlations in the WF's and consequently an acceptable convergence for the first class is reached only for Ki > 40. 
The chosen version of the Minnesota potential has the exchange parameter u = 1. As mentioned before, the point- 
Coulomb potential is included in the calculation, however, in the WF we have included only states with T = 0. In 
both cases, we observe a good agreement between the different theoretical estimates. 

The three potentials examined so far are given as functions of gaussians and thus depend on rf^ . As is well known, 
in such a case the convergence of the HH expansion as a function of the grand angular quantum number is exponential 
and fast. We actually observe such a behavior in all three cases. However, especially for the class CI, the convergence 
is relatively more difficult for the two models with a repulsive core than in the Volkov case, confirming that it is this 
class which is mostly responsible for the need of the two-body correlations. 

The next examples considered are the MT-V and MT-I/III potentials. They are given as a superposition of 
Yukawians and have a strong repulsive core with a l/r divergence. As already mentioned they represent the most 
challenging problem for the HH expansion, due to the difficulty of constructing accurate two-body correlations at 
short interparticle distances, where the cancellation between kinetic and potential energy is critical. As can be seen 
by inspecting Table EI the BE for the MT-V is slightly underestimated. We have already discussed this case in the 
previous subsection and we have seen that it is possible to obtain very precise estimates for the "missing" BE using 
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the known behavior A (x 1/ K . Adding this "missing" BE to the value B = 31.347 MeV brings the HH results 
very close to the estimates computed by other techniqueis. For the (s-wave) MT-I/III we observe that our estimate 
is already close to the very precise calculation of Ref. 18]. The "missing" BE in this case is estimated to be 21 keV, 
bringing our estimated BE to be 30.331 MeV. 

Let us now consider the calculations performed using the realistic models of the NN interactions (see Table IV1|I . 
Again the value fi^/m — 41.47108 MeV fm^, corresponding to 2/m — l/nip + l/m„, has been used. Let us consider 
first the calculations performed without any 3N interaction. We have considered here the AV18 and the Nijmegen II 
(Nijm-II) interactions models. Both potentials belong to the group of the modern NN potentials which reproduce 
the NN Nijmegen data set with a P^r datum « 1. They have been selected since they are local in coordinate 
space, while other modern potentials either have a "non-local" term like (Nijmegen I potential 0), or are given 
in momentum space (Bonn interaction 0). Note that our technique does not, in principle, present any difficulties in 
treating these other kind of potentials. The only problem is that now it is not possible to solve the hyperradial second 
order differential equations by the method proposed in Ref. . Work is in progress to overcome this difficulty and 
to compute the A = 3 and 4 WF's also with non-local potentials in coordinate or momentum space. 

The convergence of the HH expansion in the case of the AV18 potential has been already discussed in the previous 
subsection. An analogous pattern of convergence is also found for the Nijm-II potential. In Table IVII the results 
for the BE and other properties are compared with the results of other techniques. Note that in the Nijm-II model 
we have included also the electro-magnetic interactions, in addition to the Coulomb potential, as in the case of the 
AV18 potential. These terms contribute an additional —0.07 MeV to the BE and this explains the difference with the 
reported FY calculation, where they were not included. By taking into account this fact, our Nijm-II BE agrees well 
with the corresponding value obtained using the FY equations. Moreover, by taking into account the "missing" BE 
estimated as explained previously, our results practically reproduce the FY ones, by again taking into account the 
quoted 50 keV uncertainty of the latter method (sj ■ 

We now consider the inclusion of the 3N interaction. We have considered here two models: the already discussed UIX 
and Tucson-Melbourne TT] (TM) model. In the latter case, we have used the modified version TM', more consistent 
with chiral symmetry 60], with the cutoff parameter fixed to be A = 4.756 TOtt 3]. We have used them together 
with the AV18 potential (AV18-hUIX and AV18+TM' models). The cutoff parameter of the TM' 3N interaction was 
chosen to reproduce the BE of ■^He. The inclusion of the UIX or the TM' models of the 3N interaction does not 
change the convergence behavior of the HH expansion and also in these cases it is possible to obtain nearly converged 
results (they have been obtained in the approximation described in Subsec. lTVD|l . Note that in the AV18-I-UIX case, 
the HH and FY estimates for the BE are slightly above the GFMC result. This is probably due to fact that in the 
GFMC technique, the and (L • S)"^ terms of the NN interaction are not treated exactly and therefore the GFMC 
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estimates have to be regarded as an upper bound of the true ground state energy. 

In summary, the HH expansion has proved to be flexible enough to describe accurately the a-particle bound state 
using realistic NN and 3N interaction models. 



In the calculations performed with the realistic NN and NN+3N interaction we have included components with 
total isospin T — 0, 1 and 2 in the WF. The calculated percentages of the waves with T — 1 and T = 2 for the AV18 
and AV18+UIX models are reported in Table IVlTl The results obtained by the FY calculations 51] have been also 
reported. These components appear in the WF when the class C6 is included in the expansion. From Table Hill it can 
be seen that the convergence of the BE for that class is reached without difficulties including states up to if = 16. 
However, the percentage values of the T = 1 and 2 states have been found to converge substantially more slowly and 
HH functions of class C6 up to if = 32 have to be considered. The contribution to the BE of the C6 states with 
if > 16 is very small, less than 1 keV. 

As can be seen by inspecting Table IVTll the percentages of the components with T = 1 and T = 2 in the a-particle 
wave function are extremely small. For the AV18 potential, they are in good agreement with the FY estimates |5l| . 
The percentages obtained using the Nijm-II potential differ by about 40% with respect to those obtained with the 
AV18. The inclusion of the 3N interaction tends to reduce them slightly. The adopted models of 3N interaction do 
not contain any isospin mixing term. 

The knowledge of the T = 1 and 2 percentages is important for parity violating experiments of electron scattering 
on '*He, devoted to studying the admixture of strange quark ss pairs in the nucleon. Information on this quantity 
can be extracted from the measurement of the "left-right" asymmetry Aj^ji of polarized electrons on a target nucleus, 
resulting from the interference between the electromagnetic and the weak neutral current mediating the scattering 
process. The study of the asymmetry is particularly simple in case of a {J'^,T) = (0"'",0) system, since in that case 
the number of matrix elements entering this observable is small. Moreover, the use of *He as a target nucleus is also 
favored by the fact that its first excited state is at 20.1 MeV, which allows for an easy experimental control of inelastic 
processes. Indeed, there are approved experiments at the Jefi^erson Lab |6lll6^ . 

However, the extraction of the information from the experiments could be complicated by the presence of com- 
ponents with isospin T = 1 and 2in the WF of "'He. This question was analyzed in Ref. |^ and found that the 
contribution from the T = 1 isospin mixing configurations to Aj^ji was negligible. Considering the effect of the 
Coulomb potential alone, the percentage of the T = 1 component in that work was estimated to be Pt=i = 7 x 10~^. 
From the present calculation, in agreement with the study of Ref. 51], the T = 1 component results to be is 4 times 
larger and this could be of some effect on Alr- 

It is interesting to study the origin of the T = 1 and 2 isospin admixtures to the a particle WF. To this end we 
have performed a series of calculations by removing from the Hamiltonian the different terms which induce the T > 
components. Let us write 



where Hjc is the isospin-conserving part of the nuclear Hamiltonian, He the point-Coulomb interaction, iicsB the 
charge symmetry breaking nuclear interaction (namely, the operators 15-18 in AV18), iic.m. the remaining electro- 
magnetic (e.m.) interaction (finite-size effects, vacuum polarization, magnetic moment interactions, etc) and i^A the 
term originating from the proton and neutron mass difference in the kinetic energy. This latter term has not been 
included in the solution of the four body problem and its effects have been evaluated perturbatively as explained 
below. 

By approximating the Hamiltonian with Hic only, one would get no isospin admixture at all. We have then added 
the various terms one by one to iiic, and reported the results in Table IVlTll As can be seen from that table, with 
the inclusion of the Coulomb potential He, the percentage of the T — 1 state is in rough agreement (within a factor 
2) with that estimated in Ref. (g^- The percentage of the T = 2 state is very tiny in this case. When the CSB terms 
of the AV18 NN interaction are taken into account, however, the previous picture is noticeably modified and both 
components increase, in particular the T = 2 component which becomes larger than the T = 1 one. Finally, the effect 
of Hem. is rather tiny. 

To further analyze the origin of the isospin admixture components, we have repeated the approximate calculation 



C. Origin of the T > components 



H = Hic+Hc + HcSB + iic.m. + if A , 



(4.7) 



of Ref. [6^ . Let us write 



He + HcSB + iic.m. + if A = H\"^ + h\^^ + iij- 



(4.8) 
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namely as a sum of an isoscalar, isovector and isotensor term (the isoscalar part comes from the Coulomb and e.m. 

potentials). Let us treat Hjc + Hf'^ as the unperturbed Hamiltonian, and try to evaluate the T > components 
using first order perturbation theory. Namely, 

ri>0 ^ " 

where is the unperturbed ground state and '^^n n — 1,2, . . . the unperturbed excited states of Hjc + Hf \ which 
therefore have definite values of the total isospin quantum number. In particular, ^'o has T — 0, etc. The most 
important contributions to the components T 1, 2 of S'^'^^^ would come from the lowest excited states. Following 
Ref. (see also Ref. 0)' model these states as 

p ^ , (4.10) 

(vI/o|4f^T|*o)i/2 

where fiT, T = 1,2, are excitation operators of the form 

= E4(^4*)+^-0')) ' (4-11) 

ij 

^2 = E''?;(^-W^-W"(l/3)r(z)-x(j)) . (4.12) 

ij 

The operator fii generates a state {J^, T) — (0"*", 1) corresponding to a "breathing" mode where neutrons and protons 
oscillate in counter-phase. Furthermore, VI2 generates a state {J^ ,T) = (0"*", 2) with "tensor" oscillations. The T = 1, 2 
components in the *He wave function would be given by 

1^*^^)) - l^r> ^^I'l^'^l^^"^ ^ XTl^P) , (4.13) 

where e[^^ = (*^^^ \Hjc + Hf^ 1*^^)- The percentage of the T-wave is just 100|xtP- 

In the following, we applied the procedure outlined above using the AV18 potential model. ^0 is the WF computed 
with the HH expansion excluding any states belonging to the class C6. We have found 

Eq = -24.19 MeV , e[^^ = 7.05 MeV , ^ = 29.07 MeV . (4.14) 

Let us first consider the (point) Coulomb potential He, given by 



and therefore 



i<j 



Hfh - E - f . (4.17) 



The necessary matrix elements can be readily computed with the result that 

(*f='^|i7fp='^|*o)~-100keV, (*f='^|Hj^='^|*o) » -37keV . (4.18) 

and therefore Pt=i ~ 1 x 10"'^ and Pt=2 ~ 0.05 x 10^'' confirming that the T — 2 component induced by the 
Coulomb potential is smaller than the T = 1 one. In fact, the radial dependence (and strength) of Hj^^, and Hj^^, 
are the same. However, the T — 2 breathing mode has a higher excited energy and this reduces the probability of a 

(T—2,'] 

"transition" to the state \'^\ ). The values Pt=i and Pt=2 are also in rough agreement (within a factor 2) with 
the results reported in the second row of Table IVIIII It appears that in our approximate calculation, the T — 1 and 
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2 percentages are somewhat underestimated. Note that the value Pt=i « 1 x 10~^ is also in reasonable agreement 
with the estimate of Ref. |6^ . 

As seen before, this situation is reversed when the CSB part of the nuclear interaction is taken into account. The 
CSB part of the AV18 potential is explicitly given by 

HcsB=Y.{^^'K^,j)[r.{i) + T^{j)] + V^^\t,j)[T^{t)T^{j^ , (4.19) 

i<j 

where V^^\i,j) and V^'^\i,j) are functions depending on the interparticle distance r.y and on the spin operators 
of the particles i and j. The operator [Tz(i)Tz{j) — ^T,j • Tj] in Hf'^jgg induces differences between the pp and pn 
interactions (it originates mainly from the difference between the charged and neutral pion masses). These differences 
are well established and, although rather small, are of sizable value in some observables (such as the singlet np and 
pp scattering lengths). The operator [Tz{i) + t'z(j)] in Hj^l^gg induces instead differences between the pp and nn 
interaction, too. Due to the lack of precise nn data, the magnitude of this charge independence breaking term is not 
very well known, however, its strength should satisfy Hjj^g <^ ^csb- Turning to our approximate method, we have 
found that 

(*r='^l^fc5Bl*o) « -23 keV , (*r^'^|i?fctll*o) « -84 keV . (4.20) 

As expected, the matrix elements {"i/l csbI'^o) are larger than the corresponding matrix elements of Hj (jgg- 

Adding the values given in Eqs. H4.18|l and 14.20|l . we find Pt=i = 1-5 x 10"'^ and Pt=2 — 0.6 x lO^'^ as the percentage 
of the isospin components induced by the Coulomb+CSB part of the Hamiltonian. Our approximate calculation 
seems to underestimate both component percentages (in particular Pt=2), but it is in qualitative agreement with the 
results reported in Table IVnTl 

The e.m. part of the Hamiltonian produces small effects. This is also supported by our approximate calculation. 
We have found, in fact, that 

(^'f='^|iIe.m.|«'o)~6keV, (*f='^|i/e.m.|*o) «-3keV . (4.21) 

The last effect which should be taken into consideration is the mass difference between neutrons and protons. The 
associated term in the Hamiltonian is 

'^A^Ej(s!;;-i)vf <i. (4,22) 

and therefore _ffA is an isovector operator and may increase the percentage of the T = 1 state. We have not considered 
this term in the solution of the full four-body problem but we have estimated its contribution using the approximated 
procedure outlined before. We have found that 

(vl/f =i)|i/f^=i)|*o) « -9 keV . (4.23) 

The effect of i/A is therefore rather small (it would produce a change of Pt=i of about 0.1 x 10^*). 

This study therefore roughly support the result of our full calculation, namely that the charge symmetry breaking 
terms of the nuclear interaction are responsible for the relatively large T = 2 components in the a-particle wave 
function. On the contrary, the T = 1 component originates mainly from the Coulomb interaction. The other e.m. 
terms and the neutron-proton mass difference play a negligible role. 

D. Truncation studies 

In prevision of future applications of the present technique to heavier systems, we have explored the effect of 
truncating part of the NN and 3N interaction. The aim is to find a way of simplifying the Hamiltonian, still obtaining 
very precise results, with a maximum deviation of the order of 0.1% with respect to the results obtained without any 
approximation. We have explored both the effect of neglecting the NN interaction when the total angular momentum 
j of the pair is greater than a given value jm and the effect of neglecting the 3N interaction when it acts on HH states 
with grand angular quantum number greater than a given Ky[. 

Let us first discuss the case of truncation of the NN interaction. We have considered an NN interaction which 
vanishes when acting on pairs with total angular momentum j > jm, and have varied jm to see the effect on the 
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a-particle BE. We have considered the AV18 interaction, since its operatorial form allows one to compute it for 
arbitrary values of j. The results obtained can be seen in Table HXI The calculation with jm — oo means that we 
have retained the NN interaction acting in all states. By inspecting the Table, it can be seen that by taking jm > 6 
the BE and other quantities vary very little. Therefore, it seems safe to retain the NN interaction as acting only on 
states with j < 6 8. 

Let us now consider the 3N interaction. The behavior of the radial parts of the UIX or TM' 3N potential are 
rather soft at short interparticle distances. Since the large grand angular quantum number components in the WF 
are induced by the repulsive core of the potentials, this suggests that the correlations induced by the 3N interaction 
would not need such high components. Therefore, we have included in the Hamiltonian an effective 3N interaction of 
the kind 

Wkm {h j, k) = Pl^W{i,j, k)PKM (4.24) 

where W{i,j, k) is one of the standard 3N interactions and Pxm ^ projection operator which gives when it acts on 
four body HH states with a given grand angular quantum number K if K > Km- Actually, Wk-m {hh k) is an effective 
4-body interaction. We have then studied the effects on the a-particle BE by varying Km- We have considered here 
the AV18+UIX model. The results obtained can be found in Table 1x1 where for simplicity we have restricted the HH 
basis to include only the first three classes (in any case, the other classes include HH states with K < 30). As can be 
seen by inspecting the Table, the BE and the other quantities depend very little on Km- Already for Km = 20, the 
corresponding BE differs from that obtained in the non-truncated case by less than 0.1%. The calculation of the 3N 
potential matrix elements between states with if < 20 30 is noticeably simpler than in the general case and still the 
results are of acceptable precision. This should allow for HH calculations including 3N interactions also for scattering 
states and for heavier systems. This conclusion is also supported by the study of Ref. 'gE] of the incorporation of the 
3N interaction in the EIHH method. Note that the results reported in Tables UTTl and IVll have been obtained using 
the 3N interactions "truncated" as in Eq. and taking Km = 30. 



E. Asymptotic Normalization Constants 



The asymptotic normalization constants (ANC's) are properties of the bound state WF which can be related to 
experimental observables. They are interesting quantities from which useful information on the nuclear structure can 
be extracted. They also provide a test of the quality of the variational WF in the asymptotic region, as we shall see. 
This test will be particularly severe in our approach, as the description of the '^He WF in terms of the four-body HH 
functions in regions where the 1-1-3 or 2-1-2 clustering configurations are dominant will be difRcult. 

Let us concentrate first on the proton-triton ANC C^* of "^He, defined by 

*4(li,l2,«3) - Cfy^, W^-...i/2(2/3p.v) $g*(|^,^^,^3) , r,,^oo, (4.25) 

where the Jacobi vectors correspond to the permutation p = 1 (the index p will be suppressed in this section) and 
Tpt = \/2/3 is the distance between the '^H center of mass and the fourth nucleon. The function <i>^* is defined as 



L 

*?(li,€2,l3) = (5^L(ii)k(l,2,3)x44l 1 , (4.26) 

where Vt(l,2,3) = M^2^i3) is the ^H WF and X4 (U) is the spin (isospin) function of the fourth nucleon. In the 
previous equation, the spin 1/2 of ^H is coupled to the spin 1/2 of the other nucleon to give a "channel" spin 5* = 0, 1. 
The channel spin is in turn coupled to L to give a total angular momentum J = 0, therefore L = S- Due to the even 
parity of the *He state, the p — '^H clusters can only be in the state S — L = 0. In Eq. H4.25|l . W-jjj{2Pr) is the 
Whittaker function which behaves irregularly at the origin and decays exponentially for r ^ oo, while (3pt and rjpt are 
determined by 



fi„ = f~iB,-B.} , V = 1'^^^ . (4.27) 

where ~ 1.44 MeV fm, fi^/m w 41.47 MeV fm^ and and Bt are the ^He and '^H BE, respectively. Finally, the 
factor ■y/2/3pt in Eq. (|4.25l) has been introduced so that the ANC C^* be adimensional. 
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In order to calculate CS*, let us introduce the "^H — ^He overlap function as 



and the ratio 



f^,ir,,) = J d% d% d% <s(y|a - v)^^^%^^*4(li,^2,^3) , (4.28) 

cp*(v) = , , 3 ^^'^'^'^ . (4.29) 

(3/2)^ y2;^PF_^^^,i/2(2/3ptv) 

If ^4 is the exact ^He WF, for rpt oo the overlap function behaves as 

3 

fptirpt) ^ (J^ ' 72^ M^-,,.i/2(2/3ptV) > (4-30) 

and therefore Cpt{rpt) C^*, allowing for the extraction of the ANC. The WF has been determined by means of 
the pair-correlated HH (PHH) technique described in Ref . 27] and is believed to be very precise Hfi] . The dependence 
to the truncation level of the HH basis used to compute the ''He WF has been studied by computing the overlap 
function for the following three different choices of the maximum values of the grand angular quantum numbers of 
the six classes as defined in Subsec. IIV Al 



{KiM, K2M, K^M, KiM, K^M, KeAi} = {56, 32, 26, 16, 14, 24} , case a , 

{KiM, K2M, KsM, K^M, K^M, K^m} = {60, 36, 30, 20, 18, 28} , case b , (4.31) 

{KiM. K2M, KsM, K^M, K^auKsm} = {64, 40, 34, 24, 22, 32} , case c . 

The corresponding ratios Cpt are shown in Fig.|21by the open circles (case a), open squares (case b) and solid triangles 
(case c). The potential used to generate the WF is the AV18 interaction. As can be seen by inspection of the figure, all 
three functions Cpt start to deviate from the expected asymptotic constant behavior already for rpt > 5 fm, showing the 
difficulty of reproducing the cluster structure of the WF by means of the four-body HH functions. From the differences 
between the three ratios, the very slow convergence of Cpt{r) as a function of {Kim,K2m,K3m,K4MtK5m,Kqm} 
results to be evident. In particular, a detailed analysis has shown that the convergence is sensitive to the value of 
KiM- The ratio Cpt(r) obtained with the larger basis shows a slightly larger "plateau" around Vpt = 5 fm, allowing 
for a crude estimate of the ANC, Cg* « 1.7 

To obtain a greater accuracy in the extraction of the ANC we have followed another procedure ^3 ■ Assuming that 
^4 and ipt are "exact", it is not difhcult to show that the overlap function should satisfy the following differential 
equation 



l^fptir) + ^fptir) + (i?4 - Bt)fptir)+gir) = , (4.32) 



.2 „2 



where r = rpt and 



(r) 



/ d% d% ^(^1^1 - r) 



e2 

+ Wi24 + Wi34 + W234 

r 



rpt 

*4(li,l2,^3) , (4-33) 



Vij and Wijk being the NN and 3N potential, respectively. As r — > 00, the function g{r) 0, and the solution of 
Eq. (|4.32(l coincides with the Whittaker function, allowing for the extraction of the ANC via Eq. (|4.29(l . We have 
computed the function g{r) with the three different choices of the "^He WF ^'4 given in Eq. (|4.31|) and reported the 
results in Fig. ^ As can be seen, the function g(r) is peaked at r w 2 fm, goes to zero exponentially and depends 
slightly on the choice of ^'4. In fact, the selected HH bases are already large enough to accurately describe the p — ^H 
decomposition for r < 4 fm. For larger distances, probably g{r) is not computed accurately using our variational WF, 
but there g{r) becomes vanishingly small and the resulting effect on the ANC is negligible. This has been checked 
explicitly by solving Eq. (|4.32() (imposing the boundary conditions /pt(0) = fpt{oo) = 0) and computing Cpt{r) for 
the same three cases as before. The results for Cpt{r) are shown in Fig. 01 by the dotted, dashed and solid lines (the 
results of the latter two cases are practically indistinguishable). As expected, for r < 5 fm, the line goes through 
the ratio functions Cpt{r) computed directly via the overlap integral and reaches a constant value, corresponding to 
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Cg*, around r = 5 fm. The extraction of the ANC can be now achieved without any difficulty and the value found is 
Cf = 1.715. 

An analogous procedure have been repeated for the n — '^He ANC and the S- and D-wave d ~ d ANC's. We have 
used the definition 



*4 



rnh 



rnh 



oo 



S V ^Pdd "Pq (€l'€2)€3j 

rdd 

. r<dd l7Ta~ ^-ndd,5/2{'^Pddrdd) ^ddgi 7' e' \ 
rdd 



rdd 



where ^[ are the set B of the Jacobi vectors defined in Eq. H3.1|l corresponding to the permutation p = 1 and 



(4.34) 
(4.35) 
(4.36) 

(4.37) 
(4.38) 



In the latter, tph and (pd are the '^He and deuteron WF, respectively, and 



3 m 



Pnh — \l 2J^^^'^ ^ ' f^dd — Y "^(^4 - 2i?(j) 



2rn 



Vdd — 1-2 rj J 

fi^ Pdd 



(4.39) 



with Bh and Bd the ^He and deuteron BE, respectively. For the d~ d case, one can also estimate the distorted-wave 
parameter D2 defined by 



D2 = j dvdd rid fddirdd)/ j dvdd rdd Iddirdd) 



(4.40) 



where fddirdd) {X = S, D) are the S- and D-wave [d — d) — ''He overlap functions, respectively, defined in analogy to 
Eq. 1)4. 28|) . The results obtained have been reported in Table lXll together with some other theoretical and experimental 
estimates available for D2 (for a more complete list of references, see Ref. H^l)- The £'2'* parameter was determined 
in Ref. [6^ using an approximated method (a cluster model) which however seems to provide an estimate rather 
close to ours. This parameter is also in reasonable agreement with the experimental values reported in Table Ixll also 
considering the difficulty of the extraction of this quantity from the experimental data. 



V. CONCLUSIONS AND PERSPECTIVES 



We have studied the solution of the Schroedinger equation for the four-nucleon ground state using the HH function 
expansion. The main difficulty when using the HH basis is its large degeneracy, accordingly a suitable selection of 
the HH functions giving the most important contributions has to be performed. In this work, the HH functions have 
been divided into classes, depending on the number of correlated particles, the values of the orbital angular momenta, 
the total isospin quantum number, etc. For each class, the expansion has been truncated so as to obtain the required 
accuracy. We have applied this procedure in particular to the study of the ground state of the a-particle using a 
number of NN and NN-f3N interaction models. In all the cases, accurate calculations of the BE and other ground 
state properties, such as the asymptotic normalization constants, have been achieved. 

A similar procedure can be also applied for solving scattering problems. The calculation of the phase shifts and the 
various observables for n — ^H and p — '^He elastic scattering is now in progress and will be published elsewhere [t^ . 

The hyperspherical formalism is adequate for treating all kinds of modern potentials, except those containing a 
hard-core. We have considered here the AVIS and Nijmegen-II NN potentials and the UIX and TM' 3N interactions. 
The inclusion of the term present in the Nijmegen-I ,6| potential does not introduce additional difficulties. As 
an example, such a term was taken into account in ref. |73| where the PHH approach was used. Moreover, the 
HH method can be easily formulated in momentum space. It can therefore be applied also to the case of the Bonn 
potential although one additional numerical integration and the solution of an integral equation are then required. 
The application of the HH technique to the ^ = 3,4 systems with the Bonn potential is actually under way. 
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At present there are only a few other methods available for accurate calculations of the four-nucleon problem, in 
particular by taking into account a 3N force. There are two other important motivations behind this work. The first 
one is to show in details that the HH expansion applied to the four-nucleon bound and scattering problems is very 
powerful even for realistic NN interactions. The second motivation is the possibility of the extension of the method to 
larger systems. The feasibility of such an application would require the solution of the following different problems. 
First, the calculation of the generalized "Raynal-Revai" coefRcients, namely, of the coefficients relating HH functions 
constructed with different sets of Jacobi vectors. The direct generalization of the algorithm proposed in Ref. 45] is 
adequate for A = 5 8. Otherwise, different algorithms could be used |4^E3i0|- Second, the computation of the 
matrix elements of NN and 3N interactions, which can be reduced to the evaluation of low-dimensional integrals as in 
the present case. In particular, the possibility of approximating the 3N interaction as acting only on HH functions of 
low fC, as discussed in Sec. lIVDl should appreciably simplify this task. Finally, the choice of an optimal subset of HH 
functions. As A grows, the number of HH states for a given K increases very rapidly. The criteria for selecting the 
subclasses of HH functions chosen in the present paper can be readily generalized to systems with A > A. However, 
additional properties of the HH function could be exploited to further reduce the number of terms in the expansion. 
For example, one could take into account the symmetry of the space part of the states constructed as a product of 
the HH functions and the spin-isospin states. Another possibility to be explored is to include classes of HH functions 
constructed with those Jacobi vectors pertaining to different partitions of the particles. For example, in the study of 
the d + '^H — > ^He +p reaction, the use of HH functions constructed in terms of 2 + 3 and 4+1 clusterizations should 
be very useful. 
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APPENDIX 

In this Appendix, the method used for calculating the matrix elements of a local 3N interaction operator H^n 
between the antisymmetric hyperangular-spin-isospin states defined in Eq. H2.10I) will be briefly illustrated. The 
major problem to be overcome is to achieve a sufficient numerical precision, so that the differential equations for the 
functions uklst,^i{p) defined in Eq. (|2.21|l could be solved without any numerical trouble. In general, a 3N interaction 
is written as follows 

Hm= E W^N{i;j.k) , (1) 

i<j<k cyclic 

where "YlicycUc represents a cyclic sum over induces i, j and k and W3Af(?; j, k) is symmetric under the exchange of the 
particles j and k. Therefore, the problem can be reduced to the computation of the matrix element of the operator 
W3Ar(l;2,3). Once the antisymmetric hyperangular-spin-isospin states are expanded in terms of the jj states as in 
Eq. 12.22|l . one has to compute the following integrals: 

W.y{p) = Jdn (S™(1,2,3,4)| Wsn{1;2,3) |SfS'^-'-(l, 2, 3, 4)) , (2) 

where the states S^-^'^'^(l, 2, 3, 4) are defined in Eq. (|2.23() (hereafter all the Jacobi vectors are chosen to correspond 
to the permutation p = 1, and therefore the index p will be omitted). In the case of the Urbana or TM-like 3N 
interactions, W3Ar(l; 2, 3) can be taken to have the general form ff^ 

6 

I^3w(l;2,3) = ^ Fp,,(ri2,ri3,Mi2,i3) Of^ , (3) 

p.,q=l 

where /ii2,i3 — ^^12 • ^13 and Vij is the relative distance between the particles i and j. In the latter equation, 0^~^'^ 
are the operators 

Off ^ 1, (t-i ■ rj), (cr, ■ a-j), (cr, • (Tj){t, ■ t^), rf^S^j, rf^S.^{T, ■ Tj) , (4) 

where Sij is the tensor operator (the factor rfj has been included in the definition of 0^=^,6 -j-j^g^^ these operators 
are polynomials in the Cartesian coordinates of the particles). 
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Since W3]v(l;2,3) depends only on the variables p,Lp3,Lp2,£,2,£,3 we can easily integrate over the variables ^i. 
Moreover, by evaluating the spin-isospin traces and the integrals over the angles ^2, ^3 (except for fi = £,2 ■ £3) one 
reduces the matrix element given in Eq. to an integral of the type: 



(/o) = / dz + z / dx \/T~x^ d^ Fp^q{ri2,ri3,^i2^i3)'Pp^q{z,x,fi) (5) 



P,<! 



1 



where 



and 



J.2 t2 

cos2^3 = 2^-l, a; = cos2(p2 = 2-4 - 1 , (6) 



Vp,,{z,x,p) = ^= J dl (S™(1, 2,3,4)10^2 0?3|S™-(1,2,3,4)) . (7) 

In Eq. 10 the integration over the angles ^2, ■^3 (except for p — £2 ■ Cs) and the trace over the spin-isospin degrees of 
freedom is implicit. This latter part of the calculation can be performed analytically in terms of Wigner _D-matrices 
and Clebsh-Gordan coefficients. The remaining two-dimensional integration over d^i = dcos^i d0i in Eq. jT)) can be 
easily performed by taking into account that the integrand is a polynomial in cos^i and cos^i of degree K + K' . 

The functions V can be therefore calculated exactly using an appropriate Gauss integration formula with a small 
number of points. On the other hand, the functions Fp^q entering the tri-dimensional integral |(SJ) are very complicated 
functions of the variables z,x,p. Therefore, the integral 101 requires the use of extended and dense integration grids 
(about (1000)'^ points) so as to yield the needed accuracy. Since, the same integration has to be repeated for each v, 
u' , the complete calculation oiW^.u' could be very time consuming. 

However, the function Vp^q{z, x, p) can be written in general as 

'Pp,q{z,X,p) = 'Pp^q{z,X,p.) + Vl+X\/l~Z^ Vp^q^Z , X , fi) , (8) 

where Vp q and q are polynomials in z, a; and fi of maximum degree N — K + K' + 2, the 2 coming (eventually) 

from the factor rfj multiplying the tensor operators in Eq. Q). More precisely T"^ (\/l + xy/l — z^V") is the even 
(odd) part of V with respect to the variable ^i. 

Now, if p(t) is a polynomial of degree n with respect to the variable t and its value in each of n-l- 1 points ti, . . . , 
is known, using the following "Lagrange interpolation" formula, p{t) can be computed exactly for any t: 

p{t)^ y: p{t.)Li-^'\t) , L[-+'\t)^ n (9) 

i— l,...,n+l J — l,...,n-t-l * 

Therefore, once three sets of points ^i, . . . , Zn+i, Xi, . . . , Xn+i and /zi, . . . , p.N+i in the interval (—1, 1) have been 
selected, and Vp'° in the {N + 1)'^ points Zi,Xj,iik have been computed, the functions Vp^q are then known exactly 
for all possible values of {z,x, n). Since the matrix elements of the 3N interaction are needed only between HH states 
with if < 30 (and therefore max[iV] ^ 100), as discussed in Sec. lIVDl this means that in practice the functions V 
have to be evaluated only a fairly small number of times. Finally, if we evaluate 



^Itk = / dz + z J dx x"^ J d/i Fp,g(ri2,ri3,/^i2,i3) X 

xLr'\z)Lf^'\x)Lr'\l^) , (10) 
1 



"^Itk = / dz + z J dx x'^ J dfi Fp,q{r 12, ri3,fii2,i3) X 

xTTT^ v/r^T^if +i)(,)i(^+i)(,)if +i)(^) , (11) 



the required matrix elements can be obtained simply as 

AT 



p,q ij,k=l 
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where the integrals given in Eqs. (|10|l and Hll() do not depend any more on the quantum numbers v' of the HH 
states and therefore can be computed with the necessary accuracy once and for all and stored on computer disks. 
In this way, the matrix elements YVyyip) obtained via Eq. 1)12(1 are obtained very quickly (with only ^ {N + 1)^ 
operations). 
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TABLE I: Number of four-nucleon antisymmetrical hyperspherical-spin-isospin states for the case J = 0, T = and n = + 
and the selected values of the grandangular quantum number K and total angular momentum L. Mkl is the total number of 
the states defined in Eq. 12.1011 . M'j^j^ gives the number of the linearly independent states with ii + + t-z < 6. See the text 
for details. 

K L=0 L=l L=2 





Mkl 


M'kl 


Mkl 


M'kl 


Mkl 


M'kl 
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1 
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10 
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6 


1 


4 


30 


4 


45 


4 


30 


3 


6 


70 


8 


135 


12 


89 


9 


8 


140 


14 


315 


27 


205 


18 


10 


252 


24 


630 


54 


405 


36 


12 


420 


41 


1,134 


96 


721 


63 


14 


660 


59 


1,890 


160 


1,190 


102 


16 


990 


90 


2,970 


250 


1,854 


158 


18 


1,430 


128 


4,455 


375 


2,760 


236 


20 


2,002 


176 


6,435 


488 


3,960 


321 


22 


2,730 


235 


9,009 


585 


5,511 


385 


24 


3,640 


282 


12,285 


675 


7,475 


445 


30 


3,876 




9,180 




16,540 




40 


10,626 




26,565 




47,145 




50 


23,751 




61,425 




107,900 





TABLE II: Quantum numbers of the first channels considered in the expansion of the a-particle WF. See the text for details. 
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TABLE III: Convergence of a-particle binding energies (MeV) corresponding to the inclusion in the WF of the different classes 
C1-C6 in which the HH basis has been subdivided. 
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AV18 


AV18+UIX 


20 










98 Q9S 
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14.902 
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29.794 


15.992 


16.162 


40 










29.962 


16.172 


16.337 
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30 008 


16.205 


16.369 
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30.024 
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70 
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27.370 
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28.382 
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28.404 
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28.414 
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34 28 
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28.427 
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28.439 


72 
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28.444 
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34 28 
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28.461 


72 


40 


34 28 


24 


12 




24.210 


28.462 






"exact" 






31.360 


24.25 


28.50 



TABLE IV: Increments of the a-particle BE A{K), computed using Eqs. H4.5^ for the various classes i = 1, . . . , 6 and the 
MT-V and AV18 potential models. The quantities c{K,p) are defined in Eq. (14.411 and {AB)i, given by c{k,p)Ai{K), represents 
the "missing BE" for having truncated the expansion over the class i up to the given value of K — Km- Finally, the "total 
missing BE" {AB)t is computed from Eq. KEl . 









MT-V 






AV18 




i 


Km 


A,{K) [keV] 


c(if,5) 


(AB). [keV] 


A,{K) [keV] 


c{K, 7) 


(AS), [keV] 


1 


72 


0.89 


8.51 


7.57 


0.10 


5.52 


0.55 


2 


40 


0.71 


4.52 


3.21 


0.91 


2.86 


2.60 


3 


34 


0.07 


3.77 


0.26 


1.16 


2.37 


2.75 


4 


28 


0.13 


3.03 


0.39 


2.30 


1.87 


4.30 


5 


24 








1.47 


1.54 


2.26 



{AB)t 



11.43 



12.46 



24 



TABLE V: The Q-particle binding energies B (MeV), rms radii (fm) and expectation value (K) of the kinetic energy operator 
(MeV) for various central interaction models as computed by means of the HH expansion are compared with the results obtained 
by other techniques. The binding energies obtained by using the extrapolation technique described in Sect. liV'Xl are enclosed 
in parentheses. 



Interaction 


Method 


B 


(^2)1/2 


(K) 


Volkov 


HH (this work) 


30.420 


1.490 


50.319 




SVM [20] 


30.42 


1.49 






HH [35j 


30.399 






ATS3 


HH (this work) 


31.618 


1.412 


74.366 




SVM [52] 


31.616 


1.42 




Minnesota 


HH (this work) 


29.947 


1.4105 


58.086 




SVM [20] 


29.937 


1.41 






EIHH [26] 


29.96 


1.4106 




MT-V 


HH (this work) 


31.347(31.358) 


1.4081 


69.792 




SVM [22] 


31.360 


1.4087 






EIHH [26] 


31.358 


1.40851 






CRCG [23j 


31.357 








FY [17] 


31.36 








ATMS [53] 


31.36 


1.40 




MT-I/III 


HH (this work) 


30.310(30.331) 


1.4380 


66.180 




FY [18] 


30.312 







TABLE VI: The Q-particle binding energies B (MeV), the rms radii (fm), the expectation values of the kinetic energy 
operator (K) (MeV), and the P and D probabilities (%) for various realistic interaction models as computed by means of 
the HH expansion are compared with the results obtained by other techniques. The binding energies obtained by using the 
extrapolation technique described in Sect. llV'XI are enclosed in parentheses. 



Interaction 


Method 


B 


(^2)1/2 




Pp 


Pd 


AV18 


HH (this work) 


24.210(24.222) 


97.84 


1.512 


0.347 


13.74 




FY [51] 


24.25 


97.80 




0.35 


13.78 




FY [54] 


24.223 


97.77 


1.516 






Nijm II 


HH (this work) 


24.419(24.432) 


100.27 


1.504 


0.334 


13.37 




FY [51] 


24.56 


100.31 








AV18+UIX 


HH (this work) 


28.462(28.474) 


113.30 


1.428 


0.73 


16.03 




FY [51] 


28.50 


113.21 




0.75 


16.03 




GFMC [4] 


28.34(4) 


110.7(7) 


1.44 






AV18+TM' 


HH (this work) 


28.301(28.313) 


110.27 


1.435 


0.73 


15.63 




FY [51] 


28.36 


110.14 




0.75 


15.67 



TABLE VII: Percentages of the total isospin components T — 1 and 2 in the a-particle ground states for various interaction 
models. 



interaction 


method 


Pt=i [%] 


Pt=2 [%] 


AV18 


HH this work 


2.8 10"^ 


5.2 10"^ 


AV18 


FY [51] 


3 10"^ 


5 10"^ 


Nijm-II 


HH this work 


1.6 10"^ 


7.4 10"^ 


AV18+UIX 


HH this work 


2.5 lO"'' 


5.0 10"^ 



25 



TABLE VIII: Effect of the inclusion of the various isospin mixing terms in the nuclear Hamiltonian on the percentages of the 
total isospin components T — 1 and 2 in the a-particle ground states. The calculations have been performed using the AV18 
model for the nuclear Hamiltonian. For the explanation of the various terms Hic, etc, see the text. 



interaction 


Pt^i [ 


%] 


Pt=2 [ 


Vo] 


Hic 












Hic + He 


1.5 10 


~3 


0.1 10" 


-3 


Hic + He + Hcsb 


3.0 10 


-3 


4.9 10" 


-3 


Hic +Hc + Hcsb + H^m 


2.8 10 


-3 


5.2 10" 


-3 



TABLE IX: Effects of the truncation of the NN potential when acting only on pairs having total angular momentum j < Jm- 
The potential model chosen is the AVIS and the selected HH basis has {Ki,K2,K3, K4, Ks^Ka} = {64, 40, 34, 24, 0, 0} . 





B (MeV) 


r (MeV) 


Pp (%) 


Pd (%) 


4 


-24.124 


97.692 


0.344 


13.713 


6 


-24.161 


97.771 


0.345 


13.724 


8 


-24.164 


97.773 


0.345 


13.725 


10 


-24.165 


97.773 


0.345 


13.725 


20 


-24.163 


97.774 


0.345 


13.720 


00 


-24.163 


97.774 


0.345 


13.720 



TABLE X: Effects of the truncation of the 3N potential when acting only on four-body HH states of grand angular quantum 
number K < Km- The potential model chosen is the AV18-I-UIX and the selected HH basis has {Ki, K2, K3, K4, K5, K(,} = 
{64,40,34,0,0,0} . 



Km 


B (MeV) 


T (MeV) 


Pp (%) 


Pd (%) 


20 


27.351 


109.71 


0.596 


15.05 


24 


27.366 


109.67 


0.596 


15.05 


30 


27.372 


109.65 


0.596 


15.05 


34 


27.374 


109.64 


0.596 


15.05 



TABLE XI: ANC's and the parameter D2 obtained with the HH expansion and the solution of the differential equation of 
Eq. 14.321 for two potential models. The D2 parameter is defined in Eq. 14.4011 . In the third row, the theoretical estimate for 
D2'' of Ref. €8] is also shown. Finally, in the last three rows, some available experimental value for the parameter 1)2'* have 
been also reported. 



Parameter 


Cf 








Df (fm^) 


AV18 


1.72 


1.67 


1.96 


-0.209 


-0.115 


AV18+UIX 


1.75 


1.69 


1.99 


-0.277 


-0.113 


Adhikari et al. [68] 










-0.12 


Karp et al. 
Merz et al. [70] 
Weller et al. [71] 










-0.3 ±0.1 
-0.19 ±0.04 
-0.2 ±0.05 




FIG. 1: Binding energy differences for the a particle for the classes CI (circles) C2 (squares) and C3 (up triangles) as function 
of the grand angular value K (see the text for more details). The potential used is the MT-V. The curves are fitted to the large 
K part of the energy differences. 




FIG. 2: Binding energy differences for the a particles for the classes CI (circles) C2 (squares), C3 (up triangles) and C4 (down 
triangles) as function of the grand angular value K (see the text for more details). The potential used is the AV18. The curves 
are fitted to the large K part of the energy differences. 



2 




r(fm) 

FIG. 3: Ratios Cpt{rpt) as function of the p — '^H distance r^t- The ratios obtained by the direct calculation of the overlap 
defined in Eq. (14.281 with the ''He WF corresponding to the gran angular quantum numbers specified in Eq. 14. 3H are shown 
by the open circles (case a), the open squares (case b) and the solid triangles (case c). The ratios obtained by the solution 
of the differential equation defined in Eq. 14.3211 are shown by the dotted (case a), dashed (case b) and solid lines (case c), 
respectively (the dashed and solid line are almost coincident). The ''He WF were generated using the AVfS potential. 



